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ABSTRACT 

We report the discovery of a circumstellar disk around the young AO 
star, HR 4796, in thermal infrared imaging carried out at the W.M. Keck 
Observatory. By fitting a model of the emission from a flat dusty disk to an 
image at A = 20.8 jttm, we derive a disk inclination, i = 72°^ from face on, 
with the long axis of emission at PA 28° ±6°. The intensity of emission does not 
decrease with radius as expected for circumstellar disks but increases outward 
from the star, peaking near both ends of the elongated structure. We simulate 
this appearance by varying the inner radius in our model and find an inner hole 
in the disk with radius Rj n = 55±15 AU. This value corresponds to the radial 
distance of our own Kuiper belt and may suggest a source of dust in the collision 
of cometesimals. By contrast with the appearance at 20.8 /im, excess emission 
at A = 12.5 /im is faint and concentrated at the stellar position. Similar emission 
is also detected at 20.8 /im in residual subtraction of the best-fit model from the 
image. The intensity and ratio of flux densities at the two wavelengths could be 
accounted for by a tenuous dust component that is confined within a few AU of 
the star with mean temperature of a few hundred degrees K, similar to that of 
zodiacal dust in our own solar system. The morphology of dust emission from 
HR 4796 (age 10 Myr) suggests that its disk is in a transitional planet-forming 
stage, between that of massive gaseous proto-stellar disks and more tenuous 
debris disks such as the one detected around Vega. 

Subject headings: Stars: individual (HR 4796) — Stars: formation and planetary 
systems — Infrared: stars 
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1. Introduction 

It is now well over a decade since coronagraphic imaging of stellar light scattered by a 
disk around (3 Pictoris provided circumstantial evidence for the presence of an extra-solar 
planetary system (Smith & Terrile 1984). Since then, coronagraphic and adaptive-optics 
(AO) searches have continued the quest for additional examples (Smith, Fountain, & Terrile 
1992; Kalas & Jewitt 1995), but none have been successful at imaging light scattered by 
dust around similar IRAS-selected candidates (see Table VII from Backman & Paresce 
1993, for example). This failure was originally interpreted to suggest that the disk around 
f3 Pic is highly unusual (Smith et al. 1992), but subsequent identification of a few other 
main sequence A stars with fractional infrared excesses of the same magnitude has been 
used to argue that ~20% of all A-type stars pass through an early phase like that of (3 
Pic (Jura et al. 1998). If true, the failure to image other candidates is likely due to their 
greater distance; limitations in current coronagraphic techniques would prevent detection 
of even (3 Pic's disk if it were several times farther away (cf. Kalas & Jewitt 1996). An 
alternate approach is suggested by successful imaging of the thermal dust emission from the 
disk around (3 Pic at mid-infrared wavelengths (Lagage & Pantin 1994; Pantin, Lagage, & 
Artymowicz 1997). As part of an effort to take full advantage of this technique, we present 
thermal infrared imaging of the circumstellar environment of HR 4796, a more distant AO 
star (Houk 1982) thought to have a disk similar to that of (3 Pic on the basis of its IRAS 
excess (Jura 1991; Jura et al. 1993; 1995; 1998). 

At a distance of 67±3 pc (n = 14.91±0.75 mas; Hipparcos Catalog), HR 4796 is three 
and a half times farther away than f3 Pic (d = 19.3±0.2 pc; 7r = 51.87±0.51 mas), but 
exhibits a strong IRAS excess with optical depth estimated to be twice as large (Jura 
1991). An M dwarf companion is detected 7.7" (~ 500 AU) from the star and is believed 
to be physically bound on the basis of common proper motion (Jura et al. 1993). The age 
of the system is estimated to be t — 10 Myr on the basis of isochrone fitting (Stauffer et 
al. 1995; Jura et al. 1998), in keeping with its identification as an outlying member of the 
Centaurus-Lupus association (Jura et al. 1993) with average age t = 12-15 Myr (de Geus, 
De Zeeuw, & Lub 1989). This age is younger than other well-known members of the Vega 
class (t ~ 100 Myr) and suggests that HR 4796 may be somewhat transitional between 
young stars (t ~ 1 Myr) with optically thick disks (Koerner & Sargent 1995; Mannings, 
Koerner, & Sargent 1997) and "Vega-type" stars (cf. Backman & Paresce 1993). 
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2. Observations and Results 

We observed HR 4796 with JPL's mid-infrared camera MIRLIN at the F/40 bent- 
Cassegrain focus of the Keck II telescope on UT 16 March 1998. MIRLIN employs a Boeing 
128x128 pixel, high-flux Si: As BIB detector with a plate scale at Keck II of 0.14" per 
pixel. The latter was determined by scanning a star across both MIRLIN's field of view 
and that of the Keck II guide camera for which the plate scale has been measured. The 
associated field of view was 17.5". Background subtraction was carried out by chopping the 
secondary mirror at a 4 Hz rate with 8" throw in the north-south direction and by nodding 
the telescope a similar distance in the east-west direction after coadding a few hundred 
chop pairs. Observations were carried out in filters centered at A = 12.5, 20.8, and 24.5 /mi 
with widths 1.2, 1.7, and 0.8 /im and for on-source integration times 6, 17, and 6 minutes 
respectively. Small dither steps were taken between chop-nod cycles. Source images on the 
double-differenced frames were shifted and added to make the final 64 x 64 (8.3" x 8.3") 
images. Infrared standards a Sco, a Sco, and a Lyr were observed in the same way at 
similar airmasses. An additional sequence of images was obtained by alternating rapidly 
between the 12.5 and 20.8 /im filters with no telescope offsets in order to accurately register 
images at both wavelengths. 

A two-color rendition of the emission from HR 4796 at A = 12.5 and 20.8 /im is 
displayed in Fig. 1. The underlying images are shown separately in Fig. 2a and 2b. Peak 
12.5 /im emission is compact and centered on an elongated source at 20.8 /im. At A = 
12.5 /im, peak flux density of 184±18 mJy comprises nearly all of the point-like emission 
rising above a low-level plateau. Fainter emission accounts for an additional 40 mJy and is 
slightly extended to the NE and SW. We estimate the stellar photospheric emission to be 
122 mJy at this wavelength by comparison with flux densities for a Lyr at lambda = 10 /im 
(Tokunaga 1998) and assuming similar temperatures, My = 5.80 for HR 4796, and 0.03 for 
a Lyr. This implies a total excess that is nearly equal to the photospheric emission; half 
of this arises from very near the star. In the corresponding image at 20.8 /im, emission is 
elongated ~3" at PA 30° with total flux density 1.88±0.2 Jy, consistent with the 20 /im 
bolometer measurement of Jura et al. (1993) (1.86 Jy). Measurements at A = 24.5 /im look 
roughly similar and have total flux density 2.27±0.7 Jy. Surprisingly, the center of the 
20.8 /mi structure is not coincident with the peak emission but lies between bright peaks 
located at each end. Below, we interpret this morphology as arising from a nearly edge-on 
circumstellar disk with an inner hole and investigate its properties with the aid of a simple 
model. 
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3. Modeling and Discussion 

Our measured flux densities are listed in Table I and plotted in Figure 3 together with 
color-corrected IRAS measurements and an upper limit obtained with the JCMT at A = 
800 /im by Jura et al. (1993). The spectral distribution of the flux densities exhibits a 
nearly black-body shape which peaks at A ~ 60 /zm, well away from the maximum stellar 
photosphere emission. It can be reproduced by emission from a model disk with mass 
approximately 1.0 M e which is largely devoid of radiating particles inside a radius of a few 
times 10 AU (Jura et al. 1995; 1998). In order to uniquely constrain the spatial distribution 
of dust, however, it is desirable to gain the maximum amount of information from images 
which resolve the circumstellar material. 

The optically thin radiation from an annulus of width dr and radius r in a flat dusty 
circumstellar disk is estimated by Backman, Gillett, & Witteborn (1992) to be 




where r ro is the geometrical optical depth perpendicular to the disk plane at a fiducial 
radius ro, T p is the particle temperature, e\ is the particle radiative efficiency relative to 
a black body, and D the 67 pc distance to HR 4796. For a distribution of grains with 
effective size, Ao, we assume e\ = 1 for A < A and E\ = (A/Ao) -1 for A > Ao (see Appendix 
D of Backman et al. (1992) for the relation between Ao and a particle size distribution). 
The grain temperature is then T p (r) = 468(L*/A ) a2 (r/1AU)~ - 4 where is the stellar 
luminosity in solar units. HR 4796 has a temperature similar to a Lyr with My = 5.80 at 
D = 67 pc. Then, assuming L* = 54 L Q , My = 0.03, and D = 8.1 pc for a Lyr, we find 
= 18.1 L & for HR 4796. Parameters which remain to be determined include Ao, power-law 
index of the optical depth 7, inner and outer radii r in and r out , r ro , inclination angle 1 (from 
face on), and position angle 9. 

The morphology of the 20.8 /iin emission in Fig. 2b provides a useful constraint on r in , 
l, 9, and r ro under the assumption of a single power-law structure to the radial optical depth. 
However, it does little to constrain r out , 7, and Ao- The latter can be better determined 
from fits to the flux density distribution, especially if values of the former can be derived 
from a fit to the image. In both cases, we choose to evaluate the uncertainties in the range 
of acceptable parameter values by calculating the probability of the models given the data. 
We use the above formulation to simulate both the image and the flux density distribution 
and calculate the reduced x 2 over the range of parameter values. The associated probability 
of the data given the model is taken to be P = e~ x . The probability of the model given the 
data is then estimated by summing and normalizing the probabilities over the whole range 
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of parameter space considered. For the estimate of t, individual probabilities are scaled by 
cos(t) to account for the a priori expectation that the system is oriented edge on (see Lay 
et al. 1997 for a more detailed description of this Bayesian approach). 

To simulate the image of HR 4796 at 20.8 /iin, we calculated the emission from a thin 
disk as prescribed above at each point of a finely sampled grid (As = 0.015") and added 
a flux density corresponding to the stellar photosphere to the central pixel. The result 
was binned to the resolution of the observations and convolved with a normalized image 
of a standard star obtained at similar air mass and within a short time of the HR 4796 
observations. The model image was subtracted from the data to derive reduced \ 2 an d P 
in the usual way. Initial disk parameter values were estimated from the distribution of flux 
densities in order to determine orientation parameters t and 9. The most likely values of 
these were then taken in a second model fit while varying r in , r ro , and 7. Assuming the best 
values from the latter, a repeat fit to 1 and 9 was completed as a check. The final results 
for l and 9 agreed with the initially derived values 1 = 72°tgo and 9 = 28° ± 6°. 

To test for the presence and size of any clearing in the disk, we considered a range of 
values for 7 and r ro , since these also affect the radial intensity of emission. The value of 7 
was varied between and -2.5; r ro was considered over the range 10~ 5 < t 50 au < 1- The 
appearance of the 20.8 /zm emission is not sensitive to increases in the choice of outer radii 
larger than the boundary of emission in the image, since the temperature rapidly decreases 
to levels at which emission is negligible at that wavelength; we adopted r out = 125 AU. The 
resulting estimate of r« n is evident in the plot in Fig. 4. It is clear from this figure that 
the source of emission in Fig. 2b cannot extend inward to within 20 AU of the star as an 
extrapolation of a single power-law profile to the radial optical depth. In fact, the inner 
radius of the disk is estimated to be a good deal larger, r in = 55±15AU. 

The most probable model is displayed together with the residuals in Fig. 2c and 2d. 
It is apparent that, although the presence of a hole leads to a good match with the outer 
disk, some emission close to the star remains in the residual image. This emission looks 
similar to the 12.5 /xm image in appearance and may have a common origin in a warmer 
population of dust grains near the star with greatly decreased optical depth relative to that 
of the outer disk. The unresolved flux density in excess of the stellar photosphere is 62 
mJy at A = 12.5 /xm; excess residual flux density at A = 20.8 /xm is 96 mJy These values 
suggest mean temperatures of 200-300 K and associated distances from the star of 3-6 AU, 
depending on whether the grains are the same size as in the outer disk or radiate as black 
bodies. It holds in either case that the majority of grains giving rise to the emission lie well 
interior to the 55- AU radius of the hole and in a radial zone which corresponds to that of 
the zodiacal dust in our own solar system. 
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The properties of the outer disk are best estimated with measurements of the emission 
at far infrared and millimeter wavelengths which probe cold dust (cf. Beckwith et al. 1990). 
Assuming an inner radius derived from our fit to the high-resolution 20.8 /im image, we try 
to estimate these properties with a model of the spectral distribution of flux densities. In 
addition to IRAS fluxes and an upper limit at A = 800 /im, our measurements in narrow 
bands at A = 12.5, 20.8, and 24.5 jum provide new constraints on the distribution of 
material. Our image at 12.5 /im demonstrates further that much of the infrared excess at 
that wavelength does not arise from the "disk" outside r in = 55 AU. Consequently, we do 
not attempt to fit either that measurement or the IRAS 12 /im flux density. 

We varied four parameters, t 50 au, 7, r out, and A , over the ranges, 10~ 5 < t 50 au < 1, 
-4.5< 7,0, 56AU< r out <200AU, and 10/im< A < 80/im in a fit to all flux densities at 
wavelengths longer than 20 jum plotted in Fig. 3. Flux from the inner region was derived 
from a fit to the "point-source" flux densities at 12.5 and 20.8 /im and subtracted (89, 
88, 23, and 6 mJy at 24.5, 25, 60, and 100 /im). The preferred parameter values are 
r out ~ 55 — 80 AU and A ~ 10 — 40/im; t 50 au an d 7 are not well determined within these 
ranges. Unless maintained by some orbiting body, an extremely narrow ring is probably 
too short-lived to be likely. Consequently, we assume r out = 80 AU for the model plotted in 
Fig. 3, together with r in = 55 AU, r out = 80 AU, Ao = 30 /im, t^au = 0.065 and 7 = -2.0. 

4. Summary and Conclusions 

These observations provide direct evidence for the existence of a disk around the 
young AO star, HR 4796, with size and orientation that are corroborated by independent 
observations taken 2 days later (Jayawardhana et al. 1998). Further, the morphology of 
emission at A = 20.8 /im reveals an inner hole in the disk with radius 55±15 AU. Emission 
in excess above the photosphere at 12.5 /im arises predominantly from a region interior 
to this radius. When compared to excess 20.8 /im emission at this location, the 12.5 
/im radiation yields a source temperature of a few hundred K, corresponding to a radial 
distance of a few AU from the central star. The outer and inner components of the dust 
distribution are analogous, respectively, to the Kuiper Belt and zodiacal components of 
dust within our own solar system (Backman et al. 1997). Taken together, these properties 
are startlingly like those believed to have existed in the late bombardment stage of the early 
solar system, when collisions between cometesimals in the outer solar system and between 
planetesimals or asteroids in a terrestrial planet-forming zone were likely to have generated 
a large population of smaller dust grains distributed in much the same way as those around 
HR 4796. 
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Fig. 1. — Composite 2-color image of the emission from HR 4796 at A = 12.5 (cyan) and 20.8 
/im (red). Registration was accomplished by taking images at the same telescope position 
and rapidly cycling between both filters. 

Fig. 2. — a) Image of the intensity of emission from HR 4796 at A = 12.5. b) As in a, but at 
20.8 fim. c) Model of the 20.8 /zm emission from a flat dusty disk with a 55 AU inner hole 
and a central star with flux density matched to that of HR 4796. The model was generated 
as described in the text and convolved with a PSF derived from a standard star, d) Residuals 
from the subtraction of the model from the image with exaggerated intensity scale. Excess 
20.8 /im emission near the star has a flux density of 96 mJy 

Fig. 3. — Plot of the flux densities for HR 4796 listed in Table I. The solid curve marks the 
total emission from the stellar photosphere (dashed line) and a model of the emission from 
a disk with 55 AU inner hole (dotted line). 

Fig. 4. — Plot of the probability distribution for values of the radius of an inner hole in the 
disk. As described in the text, the radial power-law index of the geometric optical depth 
was also varied from 7 = -2.5 to 0.0. The distribution over all values of 7 is shown as a solid 
line. That for 7 = 0.0 is plotted as a dashed and dotted line, and for 7 = -2.5 as a dotted 
line. 
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Table 1. Flux Densities for HR 4796 



^eff 


sx 


Flux Density 


Uncertainty 


Photosphere 


Excess 


(fun) 


(H 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


12.5 


1.2 


0.223 


0.018 


0.122 


0.101±0.018 


20.8 


1.7 


1.880 


0.170 


0.047 


1.813±0.170 


24.5 


0.8 


2.270 


0.700 


0.033 


1.994±0.700 


12.0 


6.5 


0.309 


0.028 


0.136 


0.173±0.028 


25.0 


11.0 


3.280 


0.13 


0.032 


3.250±0.130 


60.0 


40.0 


8.640 


0.43 


0.006 


8.630±0.430 


100.0 


37.0 


4.300 


0.34 


0.002 


4.300± 0.340 


800.0 


100.0 


< 0.028 









■ 


\ 


* 


ft 


a) 


b) 


c) 


d) 




o 

CD 



O 



O 
CN 



LO 

d 



o 



